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i  ^AN  INVESTIGATION  INTO  THE  REAL  GAS  EFFECTS  OF _CRYOGENIC  J1ITROGEN 

IN  INVISCID _H0MENTR0PIC  FLOW, 

'  '  - - 1 

C.  M.^Albone 


SUMMARY 

As  a  contribution  to  the  investigation  of  the  suitability  of  using  cryogenic 
nitrogen  as  the  test  gas  in  a  high  Reynolds  number  transonic  wind-tunnel,  a  study 
is  made  here  of  the  real  gas  effects  of  nitrogen  at  low  temperatures.  The  study, 
which  is  limited  to  inviscid,  homentropic  flow  of  a  non-conducting  gas,  takes  the 
form  of  an  independent  confirmation  of  results  by  Kilgore  et  al^ .  A  recent  paper 
by  Wagner  and  Schmidt^  on  this  subject  employs  a  different  equation  of  state  from 
that  used  here  and  their  investigations  cover  more  than  just  homentropic  flow. 

The  new  contribution  in  this  Memorandum  is  that  the  use  of  a  simplified  equation 
of  state  enables  an  expression  for  enthalpy  (and  hence  the  terms  in  Bernoulli's 
equation)  to  be  derived  by  analytic  integration. 

A  shortened  version  of  this  Memorandum  was  presented  at  the  First  Inter¬ 
national  Symposium  on  Cryogenic  Wind-Tunnels  at  Southampton  University, 

3-5  April  1979.  _ l - - - - - - - 
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1  INTRODUCTION 

A  proposal  for  a  new  high  Reynolds  number  wind-tunnel  involves  the  use  of 
very  cold  nitrogen1  (called  cryogenic  nitrogen)  as  a  test  gas.  The  much  reduced 
viscosity  and  increased  density  at  low  temperatures  gives  such  a  facility  a  high 
Reynolds  number  capability  without  the  need  for  extreme  pressurisation  or  the  use 
of  large  scale  models.  We  wish  to  determine  to  what  extent  the  pressure  distri¬ 
bution  over  a  configuration  tested  in  cold  nitrogen  differs  from  that  for  the 

same  configuration  in  warm  air.  Some  differences  may  be  anticipated,  since  it  is 
3 

known  that  the  equation  of  state  for  a  gas  at  low  temperatures  departs  from  that 
of  a  perfect  gas. 

As  a  contribution  to  the  investigation  of  these  real  gas  effects,  a  study 
is  made  of  the  equations  governing  the  flow  of  cold  nitrogen  in  a  wind-tunnel  to 
see  how  they  differ  from  those  governing  the  flow  of  warm  air.  In  order  to 
simplify  this  study,  certain  assumptions  are  made  regarding  the  nature  of  the  flow. 
The  author  has  been  unable  to  find  any  evidence  which  shows  that  vibrational  or 
rotational  relaxation  effects  are  significant,  and  so  accordingly  these  are 
neglected,  as  too  are  the  effects  of  heat  conduction.  Since  we  wish  to  consider 
high  Reynolds  numbers,  we  will  assume  the  flow  to  be  inviscid  and  accept  that  any 
conclusions  we  may  draw  apply  to  flows  without  large  scale  separation,  and  then 
only  to  the  region  outside  the  boundary  layer.  Our  study  will  further  be  restric¬ 
ted  by  considering  the  flow  to  be  steady  and  irrotational  and  hence  (with  a 
uniform  flow  far  upstream)  homentropic*.  This  restriction  prohibits  a  consider¬ 
ation  of  flows  with  shock  waves,  although,  if  these  are  weak,  any  conclusions  may 
carry  over  in  an  approximate  way.  The  existence  of  a  velocity  potential,  however, 
reduces  the  problem  to  that  of  solving  a  second-order  partial  differential 
equation  (called  the  potential  equation)  for  the  velocity  potential  $  ,  coupled 
with  one  or  more  auxiliary  equations  for  the  speed  of  sound,  a  .  The  form  of  the 
potential  equation  is  independent  of  real  gas  effects,  since  it  is  derived  from 
the  laws  of  conservation  of  mass  and  momentum.  The  auxiliary  equations  are  how¬ 
ever  dependent  upon  real  gas  effects.  In  general,  they  consist  of  Bernoulli's 
equation  and  the  equation  of  state  together  with  a  constraint  provided  by  the 
condition  of  homentropic  flow. 


*  The  word  homentropic  is  used  in  this  Memorandum  to  describe  a  flow  in  which 
there  are  no  variations  of  entropy.  Many  readers  would  prefer  the  word  isen- 
tropic,  which  means  that  entropy  is  constant  along  a  streamline,  but  it  does 
not  imply  the  existence  of  a  velocity  potential. 


We  confine  our  investigations  to  a  study  of  these  auxiliary  equations,  and 

3 

the  starting  point  for  this  study  is  an  equation  of  state  for  nitrogen  under  a 
wide  range  of  conditions.  In  section  2,  this  equation  of  state  is  simplified  for 
the  restricted  temperature  and  pressure  range  of  a  cryogenic  wind-tunnel.  Expres¬ 
sions  for  specific  heats  and  the  speed  of  sound  are  derived  in  sections  3,  4  and  5. 
Numerical  values  for  some  of  these  parameters  under  cryogenic  conditions  are 
given  in  section  6,  together  with  an  approximate  homentropic  expansion  law  derived 
from  these.  The  exact  homentropic  flow  condition  is  derived  in  section  7,  and, 
in  section  8,  an  exact  form  of  Bernoulli's  equation  is  given.  The  near-ideal 
behaviour  of  cryogenic  nitrogen  is  illustrated  in  a  one-dimensional  expansion 
through  a  nozzle  described  in  section  9.  Some  concluding  remarks  appear  in 
section  10. 

2  A  SIMPLIFIED  FORM  OF  JACOBSEN'S3  EQUATION  OF  STATE  FOR  NITROGEN 


Jacobsen's  32-parameter  equation  of  state  for  nitrogen  is  valid  for 
pressures  in  the  range  1  atm  to  10000  atm  and  temperatures  in  the  range  65  K  to 
2000  K.  An  approximation  to  Jacobsen’s  equation  has  been  derived  (for  the 
temperature  range  100  K  ■+  300  K  and  for  the  pressure  range  1  atm  5  atm),  which 
takes  the  form  of  the  ideal  gas  equation  of  state  plus  the  first  term  of  an 
expansion  in  powers  of  p 

p  =  pRT  +  p2f (T)  ,  (1) 

1  ? 

where  f(T)  *  NjT  +  N^T 2  +  N^  +  N^/T  +  N^/T  ,  p  is  the  pressure  in  atmos¬ 
pheres,  p  is  the  density  in  moles/litre,  T  is  the  temperature  in  kelvins  and 
R  is  the  gas  constant,  taken  as  R  =  0.082.  The  coefficients  Nj-N,.  take  the 
following  values: 


N|  -  0.00136 


34.10  ; 


;  N2  »  0.10703  ; 


-  2.439 


-  4223.74  . 


At  the  extreme  conditions  of  T  =  100  K  and  p  “  5  atm,  equation  (1)  gives 
(p  -  pRT)/p  to  within  2Z  of  Jacobsen's  values  (that  is  to  say  within  2Z  of  the 
deviation  from  a  perfect  gas).  The  advantages  to  be  gained  from  having  an 
equation  of  state  in  the  simplified  form  (1)  will  be  apparent  in  later  sections 
of  this  Memorandum. 
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3  EVALUATION  OF  SPECIFIC  HEAT  AT  CONSTANT  VOLUME  (C  ) 
_ v 

4 

Equation  (38)  on  page  45  of  Howarth  implies  that 


/3C  \  „  /.2  \ 

VP  l  P2  \3T2/ 


(2) 


Using  equation  (1),  we  have 


(!ft  ■  *  >2f'  • 


(3) 


where  f'  denotes  df/dT  ,  and  so 


«  p2f" 


Thus  from  equation  (2)  we  have 


a 


=  -  Tf" 


and  so  on  integration 


Cv  =  -  Tf'p  +  g(T)  , 


where  g(T)  is  a  function  of  T  which  is  tabulated  by  Jacobsen  .  In  the  range 
T  =  100  K  -*■  300  K,  g(T)  does  not  depart  from  the  (ideal  diatomic)  value  of 
2.5R  by  more  than  0.03%.  We  accept  this  approximation  and  adopt  the  following 


expression  for  Cy: 


C  «  2.5R  -  Tpf" 

v 


4  EVALUATION  OF  SPECIFIC  HEAT  AT  CONSTANT  PRESSURE  (Cp)  AND  Y 

4 

Equation  (32)  on  page  45  of  Howarth  implies  that 


(4) 


_ _ \ 

i  nifiiiD-rTteafiYtfmha 
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C 

P 


cv  ■  -4 (HI  {& 


Since  pressure  is  a  function  of  p  and  T  ,  we  may  in  general  write 


dp  = 


Hence 


On  substituting  this  expression  into  equation  (5) ,  we  obtain 


From  equation  ( I ) ,  we  have 


RT  +  2pf  . 


Substituting  from  equations  (3)  and  (8)  into  (7),  we  have 


C  - 
P 


C 

v 


T_  (PR  +  P2f')2 
2  (RT  +  2pf)  * 


which  using  equation  (4)  becomes 


C 

P 


2.5R  -  Tpf"  + 


T(R  +  pf')2 
(RT  +  2pf ) 


The  ratio  of  specific  heats,  y  ,  may  thus  be  evaluated  as 


r 


C 

_R 

c 


using  equations  (4)  and  (9) . 
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5  EVALUATION  OF  a2 

Equation  (6)  implies  that 


where  S  is  the  entropy,  and  so  by  definition,  the  speed  of  sound,  a  ,  is 
given  by 


2 

a 


Equation  (35)  on  page  45  of  Howarth 


4 


implies  that 


(ID 


Hence 


(12) 


On  substituting  for  OT/3p)g  from  equation  (12)  and  for  (3p/3T)p  from 
equation  (7)  into  equation  (11),  we  obtain 


which  simplifies  to 


(13) 


and  using  equation  (8)  this  becomes 


a2  -  y(RT  +  2pf )  . 


(14) 
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6  A  SUGGESTION  OF  A  POSSIBLE  SIMPLE  SOLUTION  FOR  COLD  NITROGEN 
We  consider  a  set  of  values  of  T  given  by 

T  -  100  (5)  120  (10)  160  (20)  300  K 

and  a  set  of  values  of  p  for  each  value  of  T  given  by 

p  =  1  (1)  5  atm  . 

2 

For  each  value  of  p  and  T  ,  the  values  of  p  ,  y  and  a  are  evaluated 
from  equations  (1),  (10)  and  (14)  respectively,  together  with  the  compressibility 
factor  Z  defined  by 

p  =  pRTZ 


so  that 

2  -  1  +i-  (15) 

The  variation  of  Z  with  T  is  shown  in  Fig  1  for  values  of  p  =  1  atm  and 

5  atm.  In  Fig  2,  the  variation  of  y  with  T  is  shown  for  the  same  two  values 

of  p  .  Its  departure  from  the  ideal  value  of  1.4  is  considerable  at  the  lowest 
temperatures.  Shown  in  Fig '3  is  the  variation  with  T  of  a  ,  where  a  is 
defined  as 

a2  =  HE  .  (16) 

P 

For  cryogenic  nitrogen,  we  see  that  even  though  y  departs  significantly 

from  1.4  at  p  =  5  atm  as  T  is  reduced  towards  100  K,  the  quantity  a  remains 

closer  to  a  value  of  1.4.  Consider  the  consequence  of  assuming  that  a  is 

constant.  We  are  then  able  to  integrate  equation  (16),  using  the  definition  of 
2 

a  ,  to  obtain  the  homentropic  relation  as 


_E_ 

a 

P 


const. 


(17) 


Hence,  if  a  is  constant,  it  is  the  homentropic  expansion  coefficient. 
Bernoulli's  equation  for  a  steady  inviscid  non-conducting  flow  is 
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h  +  iq 


const  , 


where  h  is  the  enthalpy. 

This  would  become,  using  equations  (17)  and  (16), 

a2  2 

a  ]  +  {q  =  const  ,  (19) 

where  q  is  the  local  speed  of  the  flow. 

Hence  under  the  assumption  that  a  is  a  constant  and  that  it  takes  a  value  of 
1.4,  Bernoulli's  equation  takes  the  same  form  as  that  for  an  ideal  diatomic  gas. 
Our  assumption  has  used  numerical  results  to  suggest  a  relation  (16)  which  enables 
an  approximate  homentropic  expansion  law  (17)  to  be  postulated  resulting  in  a 
simple  form  for  Bernoulli's  equation  (19).  The  argument  for  the  existence  of  a 
simple  form  of  Bernoulli's  equation,  like  (19)  would  be  more  convincing  if 
numerical  approximations  were  not  employed  until  (if  at  all)  a  much  later  stage. 

To  this  end,  we  next  derive  the  exact  homentropic  expansion  relation  for  cryogenic 
nitrogen. 

7  THE  HOMENTROPIC  FLOW  RELATION  FOR  CRYOGENIC  NITROGEN 
For  homentropic  flow,  equation  (12)  implies  that 


Using  equations  (3)  and  (4)  for  (8p/8T)p  and  Cv  respectively,  we  have 


.  5R  -  pTf 
T 


dp  . 


which  becomes  on  integrating, 


2.5R  1 


n  T  -  j pf" 


dT  =  R  In  p  +  ff'dp 


Integrating  the  last  term  on  the  right  hand  side  by  parts  and  rearranging  gives 


!^r) 


=  0.4pf'  +  const. 
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Numerical  results  for  a  typical  homentropic  expansion  in  a  cryogenic  wind- 
tunnel  show  that  equation  (20),  when  combined  with  the  equation  of  state  (1),  is 
very  closely  approximated  by  equation  (17)  with  a  »  1.4.  This  evidence  is 
encouraging,  but  we  would  like  to  take  the  analysis  one  step  further. 

8  AN  EXACT  FORM  OF  BERNOULLI'S  EQUATION  FOR  CRYOGENIC  NITROGEN 

We  wish  to  evaluate  the  enthalpy,  h(=^(dp/p)) ,  in  Bernoulli's  equation  (18), 
for  the  equation  of  state  (1) 

p  =  pRT  +  p2f (T) 

subject  to  the  homentropic  condition  (20) 


R  =  0.4pf '  +  Kj  , 


where  Kj  is  a  constant.  Integrating  by  parts,  we  have 


h  '  It  '  [p]  '  Mp)  ‘  03  *  h dp 


Substituting  for  p  from  equation  (1),  we  obtain 


h '  03  *  If  *  *  h”  ■ 


which  rearranges  to  become 


=  +  jRTd(ln  p)  +  | fdp 


Integrating  each  term  by  parts  again  gives 


h  =  ^  +  RT  In  p  +  pfj  -  J(R  In  p  +  pf')dT 


From  equation  (20)  we  have 


In  p  +  pf 


*)dT  =  J 


2.5R  In  TdT 
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and  so  we  can  eliminate  R  In  p  +  pf'  and  obtain 


]  1 


h  =  +  RT  In  p  +  pfj  -  2.5R j In  TdT 

=  ^  +  RT  In  p  +  pf  -  2.5R(T  InT  -  T)  +  const. 

P 

Using  again  equation  (20)  and  the  equation  of  state  (1),  this  can  be  reduced  to 


h  =  3.5RT  +  2pf  -  pTf'  +  const  . 


(21) 


This  simple  analytic  form  for  enthalpy  for  cryogenic  nitrogen  does  not  appear  to 
have  been  presented  before.  We  may  now  write  Bernoulli's  equation  (18)  for  the 
homentropic  flow  of  a  gas  governed  by  the  equation  of  state  (1)  as 


3.5RT  +  2pf  -  pTf*  +  iq  =  const 


(22) 


From  equations  (4),  (9),  (10)  and  (14)  for  C  ,  C  ,  y  and  aT  respect- 
.  .  2  P  v 

ively,  we  may  write  a  as 


rt  +  5„f  +  T(R  +  Pf*)2 
RT  Pf  2.5R  -  pTf" 


(23) 


2 

Hence  we  may  introduce  a  into  Bernoulli's  equation  (with  no  approximations 
other  than  those  implied  by  the  equation  of  state  (1))  to  obtain 


h" 


const  , 


RT  +  2pf  + 


where  8  =  1  + 


T(R  +  pf')2 
2.5R  -  pTf " 


3.5RT  +  2pf  -  pTf ' 


(24) 


Values  of  8  have  been  obtained  for  the  same  range  of  T  and  p  as  had 
been  used  for  the  evaluation  of  a  in  section  6.  At  p  =  1  atm,  8  increases 
monotonically  with  T  from  1.398  at  T  =  100  K  to  1.401  at  T  =  300  K.  At 
p  *  5  atm,  8  increases  again  monotonically  from  1.387  at  T  =  100  K  to  1.403  at 
T  =  300  K.  Values  of  8  for  the  latter  case  are  plotted  in  Fig  3,  but  for 
p  *  1  atm  the  variation  is  too  small  to  be  plotted,  and  is  omitted.  It  is  noted 


12 


that  values  of  6  vary  less  than  those  of  a  in  the  cryogenic  range,  and 
nowhere  in  that  range  do  they  exceed  ±1Z  of  the  ideal  diatomic  value  of  1.4. 

9  ONE-DIMENSIONAL  FLOW  IN  A  CONVERGENT-DIVERGENT  NOZZLE 

To  illustrate  the  near-ideal  behaviour  of  cryogenic  nitrogen  in  homentropic 
flow,  we  consider  the  one- dimensional  expansion  of  nitrogen  from  typical  wind- 
tunnel  stagnation  conditions  to  a  supersonic  flow  near  to  the  saturation  boundary. 
If  A  is  the  cross-sectional  area  of  the  nozzle,  then  by  conservation  of  mass 
we  have 

pqA  =  const.  =  •  (25) 

Bernoulli's  equation  (22)  and  the  homentropic  condition  (20)  complete  the  set  of 
three  equations  for  p,  q  and  T  corresponding  to  a  given  area  distribution. 

The  elimination  of  q  between  equations  (22)  and  (25)  results  in 

/  K  \2 

3.5RT  .  2pf  -  PTf '  -  K2  , 

which  together  with  equation  (20) , 


R  In 


0.4pf '  +  K,  , 


forms  a  pair  of  coupled  implicit  algebraic  equations  for  p  and  T  as  functions 

of  A  ,  where  Kj ,  K2  and  are  constants.  These  constants  can  be  evaluated 

from  given  reference  conditions  and  then  the  equations  solved  by  a  double- 

iterative  procedure.  SI  units  are  chosen  for  this  exercise.  The  gas  constant, 

2  -2 

R  ,  becomes  296.813  m  s  ,  and  the  five  constants  in  the  equation  of  state  (1) 
become : 


Nj  =  0.17572 

N2  =  13.829 

N3  =  -  315.  14 

N.  =  4406.06 

4 

N5  =  -  545749.45. 

rhe  proposed  stagnation  conditions^  for  the  European  cryogenic  facility  were 
considered,  namely 


D„ 


T.  = 


4.4  atm 


and 


120  K  . 
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For  Che  calculation,  reference  conditions  used  to  determine  K( ,  and 
Kj  were  not  taken  as  those  at  stagnation  because  of  numerical  problems.  Refer¬ 
ence  conditions  for  p,  T,  q  and  A  were  obtained  by  an  iterative  process  which 
converged  only  when  the  target  stagnation  conditions  were  met  to  an  acceptable 
accuracy.  The  actual  values  of  Pq  and  Tq  used  in  the  calculation  are 


4.3954  atm 


119.96  K. 


The  expansion  was  carried  out  through  the  locally  sonic  conditions  to  a  Mach 
number,  M  ,  of  1.65  just  beyond  the  saturation  boundary  .  Table  1  shows  values 
of  q,  p,  p,  T  and  M  for  given  values  of  A  ,  together  with  p/Pq*  p/pQ  antl 
T/Tq  corresponding  to  A*/A  where  A*  is  the  value  of  A  at  the  sonic 
condition.  Also  shown  in  Table  1  are  the  values  of  p/Pq>  p/Pq  and  T/T^ 
corresponding  to  the  same  A* /A  for  air  as  an  ideal  gas.  In  Fig  4,  values  of 
temperature  and  pressure  extracted  from  the  table  are  plotted  against  Mach  number. 
The  departure  (D)  from  an  ideal  diatomic  gas,  defined  as 


(jl\  .  M 

vO /  nitrogen  \p0/  i 


ideal  gas  at  the 
same  area  ratio 


'[A 

VW  ni 


trogen 


is  plotted  in  Fig  5  against  Mach  number  for  pressure  and  temperature,  assuming  a 
similar  expression  for  the  departure  for  temperature.  The  departure  for  density 
is  less  than  that  for  pressure  or  temperature  over  the  whole  range. 

1°  CONCLUDING  REMARKS 

In  this  Memorandum,  it  has  been  shown  that  the  thermal  and  calorific  imper¬ 
fections  which  give  rise  to  a  significant  variation  in  y  are  such  that,  for 
the  homentropic  flow  of  a  non-conducting  gas,  Bernoulli's  equation  takes  a  form 
which  may  be  closely  approximated  by  that  for  an  ideal  diatomic  gas.  This  form 
of  Bernoulli's  equation  has  been  derived  using  a  simplified  form  of  the  equation 
of  state,  which  enables  an  expression  for  enthalpy  to  be  derived  by  analytic 
integration.  The  behaviour  of  nitrogen  under  typical  cryogenic  wind-tunnel 
conditions  is  illustrated  by  a  calculation  of  the  one-dimensional  flow  through 
a  convergent-divergent  nozzle.  The  departure  from  ideal  diatomic  gas  behaviour 
is  less  than  0.6Z  over  the  whole  range  of  interest. 
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